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Can N-body systems generate periodic gravitational 

waves? 



T. Chiba, T. Imai and H. Asada 



o 
o 

(N 

X> 
<D 
IX, 

00 

> 
m 

ON 

o 

O 



9 L,: 
6 



Faculty of Science and Technology, Hirosaki University, Hirosaki 036-8561, Japan 



Accepted Received 



X 



ABSTRACT 

None of N-body gravitating systems have been considered to emit periodic gravita- 
tional waves because of their chaotic orbits when N = 3 (or more). We employ a 
figure-eight orbit as a specific model for a 3-body system in order to illustrate that 
some of triple stars are capable of generating periodic waves. This illustration would 
imply that a certain class of N-body gravitating systems may be relevant to the grav- 
itational waves generation. We show also that the total angular momentum of this 
3-body system is not carried away by gravitational waves. A waveform generated by 
this system is volcano-shaped and thus different from that of a binary system. Finally, 
by evaluating the radiation reaction time scale, we give an order-of-magnitude esti- 
mate of merging event rates. The estimate suggests that figure-eight sources, which 
require carefully prepared initial states, may be too rare to detect. 
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1 INTRODUCTION 

Gravitational waves represent one of the great challenges 
of fundamental physics and will open a new window in 
the astronomy (e.g. Misner et al. 1973, Centrella 2003, Mio 
2006a, 2006b). No one has ever detected directly the ripples 
in a curved space-time generated by accelerated masses. A 
(quasi-)periodic source is the most promising candidate for 
the first detection. Hence, (quasi-)periodic sources have been 
studied for the last few decades in numerous articles, where 
they assume a single star in rotation/oscillation (e.g. Kokko- 
tas and Schmidt 1999, Stergioulas 2003) or double stars in 
binary motion (e.g. Asada and Futamase 1997, Blanchet 
2006). In particular, event rates of an inspiraling and finally 
merging compact binary (Curran et al. 1995, van den Heuvel 
et al. 1996, Arzoumanian et al. 1999, Kalogera et al. 2001, 
Kim et al. 2003, Burgay et al. 2003) are discussed by tak- 
ing account of large-scale detectors such as LIGO, VIRGO, 
GEO600, TAMA300 and LISA (Centrella 2003, Mio 2006a, 
2006b). Furthermore, gravitational waves generated by a bi- 
nary system have been intensively studied by perturbations 
such as the post-Newtonian approximation (e.g. Asada and 
Futamase 1997, Blanchet 2006, Faye et al. 2004, Berti et 
al. 2006), or by full numerical simulations for neutron stars 
(Shibata 1999) and for black holes (Pretorius 2005, Cam- 
panelli et al. 2006b, Baker et al. 2006, Diener et al. 2006). 
However, much less attention has been paid to N-body grav- 
itating systems because of their chaotic orbits when N > 3. 
There are existing works which are limited in the sense that 
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they study a binary system affected by the existence of a 
third body, namely perturbations induced by the third ob- 
ject (Iokaet al. 1998, Wardell 2002, Campanelli et al. 2006a). 
Here, one may ask whether or not N-body systems can gen- 
erate (quasi-)periodic gravitational waves. The purpose of 
this paper is to answer this question. 

The main results of this paper are two; (1) For the first 
time as far as we know, we show that some of triple stars 
are capable of generating periodic waves if a certain condi- 
tion is satisfied. This illustration would imply that a class 
of N-body gravitating systems may be relevant to the grav- 
itational waves generation. (2) We show also that a solution 
employed as a specific model does not radiate the angu- 
lar momentum. A well-known example in which no angular 
momentum is radiated away is two non-spinning black holes 
in a head-on collision (Smarr 1979, Price and Pullin 1994), 
where this system is axisymmetric and thus has a rotational 
Killing vector. Therefore, we can easily understand no angu- 
lar momentum loss. On the other hand, the specific system 
employed in this paper has no Killing vector as seen below. 

Our paper is organized as follows. In this paper, we re- 
strict ourselves within 3-body systems for convenience. In 
section 2, we assume the Newtonian gravity for motion of a 
3-body system. In section 3, the waveform and the loss rate 
of the energy and the angular momentum are estimated by 
using the so-called quadrupole formula for simplicity. Fi- 
nally, we discuss the the time scale for radiation reaction in 
order to estimate merging event rates. 
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PERIODIC SOLUTIONS FOR A 
THREE-BODY SYSTEM 



It is impossible to describe all the solutions to the 3-body 
problem even for the 1/r potential. The simplest periodic 
solutions for this problem was discovered by Euler (1765) 
and by Lagrange (1772). The Euler's solution is a collinear 
solution, in which the masses are collinear at every instant 
with the same ratios of their distances. The Lagrange's one is 
an equilateral triangle solution in which each mass moves in 
an ellipse in such a way that the triangle formed by the three 
bodies revolves. Built out of Keplerian ellipses, they are the 
only explicit solutions. In these solutions, each mass moves 
on an ellipse (a circle for the equal mass case). Therefore, 
the associated gravitational waveform is a superposition of 
a waveform for each mass orbiting in an ellipse. 

An interesting solution that three bodies move periodi- 
cally in a figure-eight was found firstly by Moore by numer- 
ical computations (Moore 1993). The existence of such a 
figure-eight orbit was proven by mathematicians Chenciner 
and Montgomery (Chenciner and Montgomery 2000). This 
solution is shown to be stable in the Newtonian gravity 
(Simo 2002, Galan et al. 2002). The figure-eight seems 
unique up to scaling and rotation according to all numer- 
ical investigations, and at the end its unicity has been re- 
cently proven (Montgomery 2005). Many efforts have been 
paid to reveal some properties such as the convexity in the 
figure-eight solution (Kapela and Zgliczynski 2003, Fujiwara 
and Montgomery 2005). Furthermore, it is shown numeri- 
cally that fourth, sixth or eighth order polynomial cannot 
express the figure-eight solution (Simo 2002). Nevertheless, 
no analytic expressions in closed forms for the figure-eight 
trajectory have been found up to now. Therefore, for the 
purpose of investigating gravitational waves, we numerically 
prepare the figure-eight orbit. 

For simplicity, we assume a 3-body system with each 
mass equal to m. Without loss of the generality, the or- 
bital plane is taken as the x — y plane. The position of each 
mass is denoted by (xa, da) for A — 1,2, 3. Figure 1 shows 
the figure-eight orbit, where we take the initial condition 
as {x u yi) = (-2:2,-3/2) = (0.970,-0.243) and (0:3,3/3) = 
(-2xi, -2yi) = (-2x2, -23/2) = (-0.932, -0.865) , where a 
dot denote the time derivative (Simo 2002). When one mass 
arrives at the knot (centre) of the figure-eight, I is defined as 
a half of the separation between the remaining two masses. 
It is convenient to use i instead of a distance between the 
knot and the apoapsis, because the inertial moment is ex- 
pressed simply as 2ml . Clearly this system has no Killing 
vector as seen in Fig. 1. 

The orbital period is estimated as T = 
6.33(Gm)- 1/2 l 3/2 « W 4 (M Q /m) 1/2 (£/R e ) 3/2 sec, where 
G denotes the Newtonian gravitational constant, and Mq 
and Rq are the solar mass and radius, respectively. Figure 
2 shows a configuration of the three masses after a half of 
the period. The mass labeled by 3 starts at the knot and 
sweeps all the L.H.S. of the figure-eight during the first half 
of the period. 




Figure 1. Three masses in a figure-eight at the initial time. Each 
mass is labeled by 1, 2 and 3. The initial velocity of each mass 
is denoted by an arrow. The distance between 1 and 3 is denoted 
by t. In this plot, we take 1=1. 




Figure 2. A figure-eight at / = T/2. The velocity of each mass 
is denoted by an arrow. 



3 GRAVITATIONAL WAVES 

The quadrupole moment Iij for N masses is expressed as 
Iij = X^a-i m Ax\x 3 A , where niA denotes the ^4-th mass 
at the location x\, and i and j run from 1 to 3 (a; 1 = 
x,x 2 — y,x 3 = z). The reduced quadrupole Qij is defined 
as Qij = Iij — SijI kk /3. In a wave zone, the gravitational 
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waves denoted by hij become asymptotically 
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Figure 3. Gravitational waves from three masses in a figure- 
eight orbit. The dashed curve and the solid one denote h^J and 
h^y , respectively. This plot is in contrast to that of an equal- 
mass binary, for which gravitational waves are expressed as a 
simple sine curve with one maximum and one minimum. On the 
other hand, there is a small valley between the nearest maxima 
and a small ridge between the nearest minima in Figure 3. At the 
bottom of the valley and the top of the ridge, one of the three 
masses passes the knot of the figure-eight, corresponding to Figs. 
1 and 2, respectively. 

where r is a distance from a source (e.g. Misner et al. 1973, 
Blanchet 2006). Here, TT means the transverse (hj^ri 1 — 0) 
and traceless {hj^ — 0), where n l denotes the unit vector 
of the direction of propagating gravitational waves. Figure 
3 shows a volcano-shaped waveform from the figure-eight 
solution. 

The loss rate of the energy (E) and the angular mo- 
mentum (L % ) is given by the quadrupole formula (Misner et 
al. 1973, Blanchet 2006). Because of z\ — for each mass, 
dL x /dt and dL v /dt vanish. In a case of the figure-eight so- 
lution, we can show numerically that dL z /dt vanishes. One 
may ask why the emitted waves carry away no total angu- 
lar momentum, though each body is likely to lose its orbital 
angular momentum by gravitational waves. This apparent 
paradox happens because the total angular momentum of 
this 3-body system vanishes, while the angular momentum 
of each body does not necessarily, as seen from Fig. 1. As a 
result, gravitational waves carry away the angular momen- 
tum of each body, while it is impossible for this system to 
lose the total angular momentum as the sum. For instance 
in Fig. 1, at the initial time, the particle labeled by 3 has 
no orbital angular momentum. The particle 1 has the or- 
bital angular momentum, which is carried away by gravita- 
tional waves. It is true of the particle 2 but in the opposite 
sign. Therefore, these changes in the orbital angular mo- 
mentum cancel each other. This cancellation is expressed as 
X)a=i La = 0, where La denotes the angular momentum 
of each body. No angular momentum radiation as a con- 
sequence of such a cancellation is in contrast to the case of 
the head-on collision of two non-spinning black holes (Smarr 
1979, Price and Pullin 1994). In the head-on collision case, 
each black hole moves without any orbital angular momen- 
tum with respect to the common centre of mass. In fact, the 
angular momentum of each black hole is not carried away 
(Lb — for each black hole). Clearly no total angular mo- 
mentum is carried away. 



The energy loss rate is estimated as 
As a result, the radiation reaction time scale is evaluated as 




which gives a rough estimate of the collision time. 

It is worthwhile to mention that the dependence of T, 
dE/dt and taw on m and I is the same as that on the mass 
and separation of a binary, as implied also by a dimensional 
analysis. 

Let us consider a 'compact' system of three black holes 
just before merging. To do so, we take I as a much shorter 
distance, say a dozen of the Schwarzschild radius of each 
mass. Then, the orbital period of this system is of the order 
of a few milli-seconds. The frequency of the associated grav- 
itational waves becomes a few kHz, around which the large- 
scale interferometric detectors (LIGO, VIRGO, GEO600, 
TAMA300) is most sensitive. We obtain the amplitude of 
the gravitational waves as 

where M ac h denotes the Schwarzschild radius. 

Clearly we have to take account of relativistic correc- 
tions such as the back-reaction on the evolution of the orbit. 
If the system is secularly stable against the gravitational ra- 
diation, one might see probably a shrinking (£ < 0) figure- 
eight orbit as a consequence of decreasing the total energy 
(E < 0)and keeping the vanishing total angular momentum 
(L = 0). Indeed, a chirp signal caused by the back-reaction 
at the inspiraling phase plays a crucial role in the gravita- 
tional waves detection. Furthermore, the numerical relativ- 
ity seems necessary for evolving the system in its merging 
phase after the last stable orbit. 



4 CONCLUDING REMARKS 

A figure-eight orbit for a 3-body system illustrated that 
some of triple stars are capable of generating periodic waves. 
We showed also that the total angular momentum of this 3- 
body system is not carried away by gravitational waves. A 
waveform generated by this system is volcano-shaped and 
thus different from that of a binary system. Our result will 
call our further attention to, for instance, N=4 (or more) 
cases. 

Before closing this paper, we mention the possibility 
of gravitational waves from 3-body systems in a 'periodic' 
orbit such as a figure-eight. Clearly the possibility seems ex- 
tremely low though it is quite uncertain. A quick order-of- 
magnitude estimate is made as follows. As a new outcome of 
binary-binary scattering, the figure-eight orbit was discussed 
for presenting a way of detecting such an orbit in numerical 
computations (Heggie 2000). According to the numerical re- 
sult, the probability of the formation of figure-eight orbits 
is a tiny fraction of one percent. Here, we assume rather 
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optimistically one per galaxy, which implies 10 6 figure-eight 
systems within lOOMpc. If I distributes randomly with its 
mean of O(IAU), we use Eq. (|3]) to obtain taw ~ 10 s Gyr 
for I ~ 1AU. As a result, merging event rates are of the or- 
der of 10 -2 event/Gyr within 100 Mpc, which is extremely 
low. Therefore, the above estimate suggests that figure-eight 
sources may be less relevant to LIGO and other detectors. 
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